QED treatment of electron correlation in Li-like ions 
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A systematic QED treatment of electron correlation is presented for ions along the lithium isoelectronic se- 
quence. We start with the zeroth-order approximation that accounts for a part of the electron-electron interaction 
by a local model potential introduced into the Dirac equation. The residual electron correlation is treated by per- 
turbation theory. Rigorous QED evaluation is presented for the first two terms of the perturbative expansion; the 
third-order contribution is calculated within the many-body perturbation theory. We report accurate numerical 
results for the electronic-structure part of the ionization potential for the n — 2 states of all Li-like ions up to 
uranium. 
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o ■ 

^ ■ 



> 
o 

m 

o ■ 

: 

o . 

o : 

>>: 
43 : 

Oh- 



X 



I. INTRODUCTION 

Li-like ions are among the most fundamental many-electron 
systems. Excellent experimental accuracy achieved for these 
ions and their relative simplicity make them a very attrac- 
tive object for an ab initio theoretical description. Numer- 
ous measurements were performed for the Lamb shift for ions 
along the lithium isoelectronic sequence; some recent ones 
were presented in Refs. d IIS H HH, 0] • Large theoret- 
ical effort was invested during the last decade into rigorous 
calculations of Q ED effects for the 2pj-2s transition energies 

, which resulted in a sub- 
stantial improvement of the description of these transitions. 
However, the accuracy achieved in theoretical investigations 
remained in most cases lower than that of the best experi- 
mental results. Its further improvement is important since a 
comparison of theoretical and experimental results on a bet- 
ter level of accuracy will provide a test of QED effects up to 
second order in the fine structure constant. 

In the present investigation, we restrict ourself to an ab 
initio description of the electronic-structure part of the ener- 
gies of the n — 2 states of Li-like ions. By the electronic- 
structure part we understand the contributions that are induced 
by Feynman diagrams involving only the electron-electron in- 
teraction (thus excluding the diagrams of the self-energy and 
vacuum-polarization type) taken in the limit of the infinitely 
heavy nucleus. The electronic-structure effects correspond 
to a well-defined and gauge-invariant part of the energy and 
thus can be addressed separately. In order to obtain accu- 
rate theoretical predictions for the energy levels, other effects 
should be added to the electronic-structure part, which com- 
prise the one-loop self-energy and vacuum-polarization (see, 
e.g., the review jl8ll ), the screening of one-loop self-energy 
and vacuum-polarization j^, [Hfl, the two-loop QED ef- 
fects B17I1 . and the relativistic recoil effect [ 19, 20]. 

Within QED, a single electron-electron interaction is rep- 
resented as an exchange of a virtual photon and is given (in 
relativistic units h = c = 1) by the operator 



where D^ v is the photon propagator, whose expression in the 
Feynman gauge reads 



£V(w,Xi 2 ) = g^u 



exp (iy/ 'w 2 



(2) 



where Xx2 — |xi2 1 = |xj — x 2 | and a M = (1, ct) are the 
Dirac matrices. The electronic-structure part of the energy 
arises through exchanges of an arbitrary number of photons. 

The traditional approach to the treatment of the electron 
correlation (beyond the lowest-order case of the one-photon 
exchange, which is rather simple) is based on a simplified 
form of the interaction obtained within the Breit approxima- 
tion. It is achieved by taking the small-cj limit of the operator 
I(oj) while working in the Coulomb gauge (i.e., neglecting the 
retardation of the photon). We will indicate the Breit approxi- 
mation employed for the operator I(ui) simply by omitting the 
energy argument, I(u>) — ► /. The corresponding expression 
consists of two parts, referred to as the Coulomb and the Breit 
interaction, 
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where x = x/|x| and a = e 2 /(4tt) is the fine-structure con- 
stant. 

The basic traditional methods for the relativistic structure 
calculations are the many-body perturbation theory (MBPT) 
applied to the lithium isoelectronic sequence by Johnson et 
al. ll2lll . the multiconfigurational Dirac-Fock method ap- 
plied by Indelicato and Desclaux j22tl . and the configuration- 
interaction method applied to Li-like ions by Chen et al. l23ll . 
All these methods treat the one-photon exchange correction 
exactly and the higher-order electron correlation, within the 
Breit approximation only. They are thus incomplete to the 
order a 2 (Za) 3 m, which is the leading order for the two- 
electron QED effects. Rigorous QED calculations of the two- 
photon exchange correction were performed for n = 2 states 
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of Li-like ions in Refs. 

cam mm 

[l5ll . An exchange of 
three and more more virtual photons can presently be treated 
within the Breit approximation only, which leads to an incom- 
pleteness at the nominal order a 3 (Za) 2 m. 

The three-photon exchange correction can be either calcu- 
lated directly, as a quantum mechanical third-order perturba- 
tion correction with the same starting potential as in QED 
calculations, or inferred from the total energies obtained in 
relativistic-structure calculations, by subtracting the zeroth- 
, first-, and second-order perturbation terms evaluated sepa- 
rately. Care should be taken in the latter case since the sub- 
traction procedure involves large numerical cancellations. 

Calculations of the three-photon exchange correction for 
Li-like ions were performed previously in Refs. lfl4ll24ll . They 
were carried out on the Coulomb wave functions and the re- 
sults can be directly added to the existing values for the two- 
electron QED effects. Another approach to this problem is ad- 
vocated by Sapirstein and Cheng lfl3ll and consists in starting 
the perturbative expansion with a local model potential, which 
incorporates a part of the electron-electron interaction effects. 
By a proper choice of the model potential, one can signifi- 
cantly accelerate the convergence of the perturbative expan- 
sion. In Ref. fl3ll . this approach was applied to the case of Li- 
like bismuth and the results are in good agreement with those 
obtained with the Coulomb wave functions Ull [l2ll . Advan- 
tages of this approach should become more evident for the 
smaller values of Z, where the convergence of the perturba- 
tive expansion becomes slower. 

In the present investigation, we adopt the method of 
Ref. lfl3ll and apply it to the study of the electronic-structure 
effects for ions along the isoelectronic sequence of lithium. 
The approach can be regarded as a successor of the MBPT 
treatment applied to Li-like ions by Johnson et al. 12111 : the 
difference is that we perform a rigorous QED evaluation of 
the second-order correction and employ a different potential 
for the zeroth-order approximation. 

The paper is organized as follows. The choice of different 
model potentials to be used for the zeroth-order approxima- 
tion are discussed in the next section. In Sec. [Till the per- 
turbative expansion for the electron-correlation effects is con- 
structed within the MBPT approximation. We present results 
of numerical calculations for different model potentials in the 
case of lithium and compare the convergence of the result- 
ing expansions. The rigorous QED calculation of the electron 
correlation through the second order of perturbation theory is 
presented in Sec. [IV] We combine the QED values for the 
two-photon exchange correction with the MBPT results from 
the previous section to obtain the electronic-structure part of 
the ionization potential for the n = 2 states of atoms along 
the lithium isoelectronic sequence. In the last section, we col- 
lect all theoretical contributions available for the 2pi/ 2 -2s and 
2p 3 / 2 -2s transition energies of several Li-like ions, compare 
them with experimental results, and discuss perspectives for 
further improvement of theoretical predictions. 



H. CHOICE OF POTENTIAL 

We consider here three local model potentials that are sup- 
posed to account for a part of the interaction between the va- 
lence electron and the closed core. The simplest choice is the 
core-Hartree (CH) potential defined as 



f°° 1 

V C n(r) = V nnc (r) + a / dr' — p s c c (r') , 
Jo r > 



(6) 



where V nuc (r) denotes the nuclear potential (i.e., the Coulomb 
potential induced by an extended nuclear-charge distribution), 
r> = max(r, r'), p c is the density of the core electrons, 



Pc(r)=y>j c + 1) [G 2 c (r)+F?(r)} 



(7) 



the superscript "s.c." indicates that the density has to be calcu- 
lated self-consistently, j c and n c are the angular-moment and 
the principal quantum number of the core electrons, and G c 
and F c are the upper and the lower radial components of the 
wave function. The CH potential plays a special role for alka- 
line ions since a well-defined part of the screening effects can 
be accounted for by perturbing the first-order QED corrections 
with this potential (see, e.g., Ref. |8fl). Owing to this property, 
the CH potential was frequently employed for an approximate 
treatment of the screening of QED effects HQ El; a similar 
potential (without self-consistency) was used in Refs. j26ll27ll . 

The second model potential is based on results of the 
density-functional theory (DFT) and referred to as the Kohn- 
Sham (KS) potential. The local potential derived from DFT is 
given by i28ll29l[30ll 
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where pt(r) = p c {r) + Pv{r) is the total charge density, 
p v (r) = G 2 (r) + F 2 (r) is the charge density of the valence 
electron, and x a E [0, 1] is a parameter. This potential has 
a non-physical limit at large r and need, therefore, to be cor- 
rected Bill to yield 
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where N c the number of the core electrons. In our calcula- 
tions, we smoothly restore the correct asymptotic behavior of 
the potential by adding a damping exponent. The result is 
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where the parameter A has to be sufficiently small in order not 
to change the original DFT potential significantly in the region 
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of interest, but sufficiently large in order to restore the proper 
asymptotic at the cavity radius. (In actual calculations, we put 
our system in a cavity with the typical radius R = 80 /Z a.u.) 
The following value for the parameter A was employed, 
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where 7 = \J n 2 — (Za) 2 and k v and n r are the Dirac 
quantum number and the radial quantum number of the va- 
lence state v, respectively. The parameter x a was set to be 
x a = 2/3 I28H30I1 . if not specified otherwise. 

The third model potential is obtained from the numerical 
solutions of the standard Dirac-Fock (DF) problem; it is re- 
ferred to as the local Dirac-Fock (LDF) potential. The scheme 
for construction of this potential was developed previously in 
Refs. 1321 13311 . For completeness, we reproduce it here, indi- 
cating some modifications introduced in this work. 

The solution of the radial DF equations for the valence state 
v (achieved in this work with help of the GRASP package 
BMP provides us with the eigenvalue e^ F and the upper and 
the lower radial components denoted by G® F and F® F , re- 
spectively. Let us now imagine that they are the eigenvalue 
and the solutions of the radial Dirac equations with a local 
potential V v , 



d n, 



" 1 pDF 



dr r J 

= (^ F -™)CW- (12) 

(l + T) GT(r) + V v (r)F^(r) 

= ( £ ° F + m)F° F (r). (13) 



Generally speaking, one cannot find a local potential V v ex- 
actly satisfying both radial equations. We can try, however, 
to solve this problem approximately. For the states v with the 
wave function without nodes [in our case, such are the 2p 1 / 2 
and 2p 3 / 2 states], we obtain our local potential by multiply- 
ing Eqs. ( fl~2l ) and ( fl~3l by G° F and F^ F , respectively, adding 
them together, and then inverting the equation with respect to 
V v . The resulting potential is 
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where v is the valence state with n r = 0. In the case when 
the valence wave function has nodes (n r > 0), the potential 
(TPfl i has to be modified since the density p^ F (although posi- 
tively defined everywhere) vanishes in the nonrelativistic limit 
at the nodes of the upper component, making the potential to 
be nearly singular at these points. 

In order to construct a smooth version of the potential ( fl~4] i 
for n r > 0, we define the density ~p v that has an admixture 



of the density of core states with the same value of k as the 
valence state, 



DF 
w K,n Pk.ti 1 



n<n v , k=k v 



1. (15) 



By a proper choice of the weights w K . n , the density ~p v can be 
made to behave more smoothly than p v . At the same time, the 
admixture of the core density should be small in order not to 
disturb the potential outside of the vicinity of the nodes. We 
achieve this by choosing the weights to be dependent on r, 



]Texp [-10(Za) 2 (r-n) 2 )] 
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(16) 

for n < n v , and w K „,n„ oc 1, where r, are the nodes of the 
upper component. The choice of weights w K _ n in the form 
( fl~6b is the only difference of the present LDF potential from 
the one constructed in Refs. [32, 33]. 

In order to build a potential with the density ~p v , we first 
define a new function X v , representing the potential ( TBb in 
the form 



V&, v {r) = V nuc (r)+a I dr'-p» F {r') 
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X v is a smooth function also for n r > 0. Analogously, we 
introduce the functions X K .„ for each core state in Eq. ( fl~5l ). 
Finally, we define our local potential as 

f°° 1 

V LDF , v (r) = V nuc (r)+a dr'—p» F (r') 
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Obviously, Vldf.u = V^ F v for n 
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It is of interest to compare the radial dependence of differ- 
ent screening potentials defined above. In Fig. Q] we present 
such comparison for the (Is) 2 2s and (Is) 2 2pi/ 2 states of 
lithium, expressing the potentials in terms the effective nu- 
clear charge 



Z eS (r) = -[V e s(r)-V nnc (r 
a 



(19) 



where V c g(r) is one of the LDF, KS, and CH potentials. Fur- 
ther discussion of different potentials will be postponed until 
the next section, where we will be able to compare the con- 
vergence of the perturbative expansions based on them. 



III. MBPT TREATMENT 

The MBPT picture can be obtained from the full QED the- 
ory by applying the Breit approximation to the operator of the 
electron-electron interaction I{uS) — > I = J Coul + jBreit^ as 
indicated by Eqs. (O-®. Using the Breit approximation in 
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FIG. 1: Radial dependence of the effective nuclear charge Z c g(r) of different screening potentials for the (Is) 2s and (Is) 2p!/ 2 states of 
lithium. 



the second and higher orders of perturbation theory, it should 
be taken into account that a summation over the negative- 
energy part of the Dirac spectrum may lead to large spuri- 
ous effects I35fl. Consistent treatment of the negative-energy 
states in such situations can be performed only within QED 
and the common approach is to restrict the summations to the 
positive-energy part of the Dirac spectrum only. 

General formulas for the MBPT corrections to the energy 
for systems with a single electron outside a closed core are 
known and can be found in Ref. ll36ll . It is our intention, how- 
ever, to obtain a different representation for them in order to 
make explicit connections with the corresponding formulas in 
full QED. 

As the zeroth-order approximation, we take the eigenval- 
ues and the solutions of the Dirac equation with the effective 
potential V c s(r) = V nuc (r) + U(r), where U(r) is the model 
potential that accounts for a part of the electron-electron inter- 
action. Since we intend to incorporate the MBPT results into 
ab initio QED calculations, the model potential U should be a 
local one, thus excluding the possibility to use the Dirac-Fock 
potential, which is the standard choice in traditional many- 
body calculations. 

In the following, we assume the electron configuration to 
be of the form (Is) 2 v, where v is the valence electron and 
Is is a core electron (also denoted by c in formulas below). 
The contributions due to the core-core interaction will be ex- 
cluded from consideration as they do not influence transition 
energies; corrections obtained in this way can be identified as 
contributions to the ionization energy of the valence electron. 

The zeroth-order energy is defined as 
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e v — m . 
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where e v is the Dirac eigenvalue corresponding to the valence 
state v. Corrections to the lowest-order energy are treated by 



perturbation theory. Since the model potential U is included 
into the zeroth-order approximation, we have to account for 
the "residual" interaction (— U) in each order of perturba- 
tion theory and the perturbative expansion goes both in the 
electron-electron interaction and in (-U) (see Ref. ll37tl for 
details). 

We will denote corrections to the energy 
where the index i indocates the order of the correction with 
respect to the electron-electron interaction, j denotes the or- 
der in U, and k = i+j is the total order of perturbation theory. 
The following shorthand notations will be used throughout the 
paper: I a bcd = {ab\I\cd} , U a b = (a, | 1 6) . By the symbols 
SW we will denote the ith-order (in U) perturbation of the 
wave function. The first two corrections are 
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where the prime on a sum denotes that the terms with the van- 
ishing denominator should be omitted. We found it important 
to keep the summation over the complete Dirac spectrum in 
the above corrections to the wave functions, not excluding the 
negative-energy states as is customary in MBPT calculations. 
Similar conclusion was previously drawn in Ref. [ 38], where 
it was demonstrated that inclusion of some negative-energy 
states drastically reduces the potential dependence of MBPT 
results for He-like ions. 

For the MBPT corrections that do not depend on the model 
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potential, we obtain 
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where ^ c denotes the momentum projection of the core elec- 
tron, P and Q are the permutation operators [(PvPc) = (vc), 
(cv)], and the sign "(+)" on a sum indicates that the summa- 
tion is performed over the positive-energy part of the Dirac 
spectrum only. For the three-electron corrections, the num- 
bers 1, 2, and 3 numerate the electrons in the configuration (in 
arbitrary order). The operator Si acts on energy denominators 
Ai, A 2 and is defined as 
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, if Ai =0,A 27 ^0, 
if Ai = , A 2 = . 



The corrections of the form AE^' 1 ^ are just the ith-order 



[in (-U)] perturbations of the Dirac energy, 

AE^=-U VV , 

AE^ = U sWvv , 



(27) 



(28) 



AE (o,m = -U 6(1)vSWv + U vv (8^v\S^v) . (29) 

The corrections of the form A-E^ppp 1 ^ are the ith- 
order perturbations of the one-photon exchange contribution 
AE^. They are given by 
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( 2 1 131 ( 2 1 2) 

Finally, A£^[ BPT is the first-order correction to the two-photon exchange contribution A£^bpt' It consists of 3 parts that 
arise as perturbations of the wave functions ("wf ), binding energies ("en"), and propagators ("pr"), 
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The corresponding corrections are given by 
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where the operator S 2 acts on energy denominators Ai, A2 as following: 
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Collecting together all corrections to a given order of perturbation theory, we obtain 
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The expressions for the second- and third-order corrections presented above contain products of the operators / and, therefore, 
include the Breit interaction to the second and even to the third order. Keeping the Breit interaction to the second order produces 
contributions of the same order as QED effects omitted [for the two-photon exchange, the order is a 2 (Za) 3 ] and thus cannot be 
regarded as ultimately wrong. However, this contribution can be shown to originate from the region of virtual excitation energies 
of order of the electron mass, where the Breit approximation is no longer valid. We thus consider it to be more correct to treat 
the Breit interaction as a first-order perturbation only. In our actual calculations, we expand each / into a sum of the Coulomb 
and Breit parts and keep terms with the Breit interaction up to the first order only. 

The MBPT results (|37Ti-(|39]l can easily be improved by accounting for the one-photon exchange correction rigorously, using 
the exact expression for the operator of the electron-electron interaction (fTJ. In this case, approximate formulas ((23), d30t . and 
(l3TT l should be replaced by the exact ones 
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where I a bcd{A) = (ab|J(A)|cd) and 1(A) is the exact op- 
erator of the electron-electron interaction given by Eq. (Q]). 

Other notations are: A a (, = e a — St, A« = U vv — U cc , 
A$ = tf 4(1) „„ -f/, (1)cc ,and/'(A) = d/(de)I(e)\ £=A . 
The resulting set of corrections to E^ ' is given by 



AE 
AE 



AE^ 

(2) 



A^-fW + A^ ' 1 ! 1 ), 



(43) 



MBPT 

(3) 

MBPT 



(2,0|2) 



uuii| tAE^+A^ 2 ! 2 ),^) 

P (3,0|3) A P (2,l|3) 
MBPT ' MBPT 

A£; (1 ' 2|3) + AS (0 ' 3|3) . (45) 



AE 



Now we have the exact expression for the first-order correc- 
tion AE^ 1 ), whereas the second-order correction still contains 
the two-photon exchange part calculated within the MBPT ap- 
proximation only. The rigorous QED treatment of the cor- 
rection AE^ 2 ' ' 2 ) will be presented in the next section. The 
three-photon exchange contribution is more difficult to eval- 
uate rigorously; presently it can be addressed to only within 
MBPT or other methods based on the Breit approximation. 

We now discuss some details of our numerical evaluation 
of the MBPT corrections. The complete set of Dirac eigen- 
states for a given local potential V c s(r) is generated by means 
of the dual-kinetic -balance basis-set method [39] with basis 
functions constructed from B splines. The computationally 
most intensive part of the calculation is the third-order cor- 
rection AE^gp 3 ^, particularly the first term in Eq. (l25l l that 
involves a quadruple sum over intermediate states. The sum- 
mation over magnetic substates was performed numerically 
for each coupling by summing combinations of the Clebsch- 
Gordan coefficients. It is possible to derive closed expressions 
for such sums (as was done for the third-order Coulomb ex- 
change in Ref. 112 111 ). However, since it does not significantly 
influence the performance of our code, we prefer to do this 
summation numerically. 

The main problem with the evaluation of the first term in 
Eq. (1231 1 is that it involves two independent (and unrestricted) 
partial-wave summations over the Dirac angular momentum 
of intermediate states, which we denote as Ki and Ki. Two 
other summations (over K3 and K4) become finite after taking 
into account selection rules of the Racah algebra. Our ap- 
proach to the evaluation of the double partial-wave expansion 
over K\ and K2 is similar to the one used in the previous cal- 
culation of the two-loop self-energy correction IHOl l4lll . First, 
we construct a matrix of elements with fixed values of \n\ | and 

M: ^|ki|,|« 2 |> \ k i\ = 1) • • ■ and \ K ?\ = 1, . . .. Then, we ex- 
trapolate to infinity the partial sums along the diagonals of the 
matrix Jfi rei u K2 i, taking typically of about 5 terms for each 



diagonal. Finally, we extrapolate to infinity the sequence of 
the sums of zth subdiagonals of the matrix, taking typically of 
about 7 terms. Significant numerical cancellations should be 
mentioned that arise in our numerical evaluation of the third- 
order correction in the low-Z region. This leads to necessity 
to employ a sufficiently large basis set in order to achieve a 
proper convergence. In actual calculations, we used a basis 
set with the number of positive-energy states varying from 50 
in the low-Z region to 40 in the high-Z region. 

An important question that arises in actual calculations is 
how one should interpret the condition of the vanishing of the 
denominator, which is used in order to exclude terms from the 
summation in Eqs. (l24l etc. and to define the case in Eqs. d26l ) 
and ( |36*l l. The point is that, for v = 2pj, an intermediate 
state with the same orbital momentum I, namely 2pj>, appears 
in the Dirac spectrum, whose energy is separated from e v by 
relativistic effects only, 6 — £2p 3/2 — e 2p 1/0 ^ (Za) 2 . As a 
result, very small energy denominators of order 5 2 may appear 
in the expressions for the third-order energy, which would lead 
to huge numerical cancellations in the low-Z region. 

The appearance of small denominators can be avoided if 
we take into account that there could be no contributions to 
the energy of the state (ls) 2 2pj arising from the intermediate 
3-electron state (ls) 2 2pj/ with a different total momentum 
J'. This means that the total contribution of such intermedi- 
ate states is zero, although the corresponding terms can arise 
at intermediate stages of the calculation. We have verified this 
statement numerically, which served as an additional cross- 
check of our numerical routine. In the actual calculations, 
we define the energy denominator as vanishing if its absolute 
value is smaller than the fine-structure difference S. This re- 
definition greatly facilitates the numerical evaluation of the 
third-order correction in the low-Z region. 

Let us now address the following questions: (i) what is 
the best local potential to use for the perturbative expansion 
and (ii) how to estimate the residual electron-correlation ef- 
fects? In order to answer these questions, we turn to the most 
difficult Li-like atom to describe within MBPT, to neutral 
lithium. For higher-Z systems, the relative weight of the elec- 
tron correlation with respect to the binding energy gradually 
decreases, so one should expect that the perturbation expan- 
sion converges faster. If we use the pure Coulomb potential 
in the zeroth-order approximation, the perturbative-expansion 
parameter is just 1/Z. It is no longer so for other potentials, 
but the general tendency remains. Obviously, MBPT is not 
the best approach to use for very low-Z systems. Extremely 
accurate results are nowadays available for lithium; they can 
be used as a benchmark to estimate the residual correlation 
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TABLE I: MBPT contributions to the ionization potential of n — 2 
states of lithium for different model potentials, in a.u. 





Coul 


CH 


KS 


LDF 


2si/2 state: 










Dirac: 


-1.12517 


-0.18310 


-0.24861 


-0.19601 


1-phot: 


1.19370 


-0.02159 


0.05362 


-0.00031 


2-ph(MBPT): 


-0.25067 


0.01158 


-0.00334 


-0.00184 


3-ph(MBPT): 


-0.00838 


-0.00939 


0.00017 


0.00001 


Sum: 


-0.19051 


-0.20249 


-0.19816 


-0.19815 


Exact: 




-0.198 159 72 




2pi/2 state: 










Dirac: 


-1.12517 


-0.12704 


-0.18537 


-0.12867 


1-phot: 


1.40602 


-0.00401 


0.06161 


-0.00109 


2-ph(MBPT): 


-0.37129 


0.00172 


-0.00734 


-0.00045 


3-ph(MBPT): 


-0.02608 


-0.00186 


0.00120 


-0.00009 


Sum: 


-0.11651 


-0.13119 


-0.12990 


-0.13029 


Exact: 




-0.130 242 69 




2p 3/2 state: 










Dirac: 


-1.12503 


-0.12704 


-0.18536 


-0.12867 


1-phot: 


1.40571 


-0.00401 


0.06160 


-0.00108 


2-ph(MBPT): 


-0.37105 


0.00172 


-0.00734 


-0.00046 


3-ph(MBPT): 


-0.02617 


-0.00187 


0.00119 


-0.00009 


Sum: 


-0.11654 


-0.13120 


-0.12991 


-0.13029 


Exact: 




-0.130 241 17 
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FIG. 2: The deviation of other calculational results from our MBPT 
values for the ground state of light Li-like ions. Squares, closed dots, 
and open dots indicate results by Chung 1 45, 46], Johnson et al. 12 ill , 
and Blundell et al. 14411 , respectively. The gray filled area denotes the 
estimation of the residual electron correlation in our calculation (see 
text). 



effects in our calculations. The best results for lithium are ob- 
tained within the method that starts with a variational solution 
of the three-body Schrodinger problem with a fully correlated 
basis set and then adds relativistic and QED corrections term 
by term as obtained within the Za expansion l42l l43ll . 

Our numerical results for the ionization energy of the n = 2 
states of lithium are presented in Table U for four potentials: 
the pure Coulomb (Coul), the core-Hartree, the Kohn-Sham, 
and the local Dirac -Fock potentials. The entry labeled "Dirac" 
contains the zeroth-order energies E^; the next 3 lines con- 
tain the corrections given by Eqs. d43l-d45ll, respectively. The 
"exact" results represent the sum of the nonrelativistic energy 
and the Breit correction in the infinitely-heavy-nucleus limit; 
they are taken from Ref . 114211 . 

It is not at all unexpected that for lithium, the perturba- 
tion theory starting with the Coulomb potential yields a very 
slowly converging series and that the results do not reproduce 
well the "exact" relativistic energies. It is perhaps more sur- 
prising that the CH potential, which gives a reasonable value 
for the zeroth-order energy, does not nevertheless yield a well- 
converging perturbative expansion. We have to conclude that 
the residual interaction is not small enough in this case. Al- 
though the CH perturbative expansion becomes converging 
for higher values of Z, this choice of the model potential does 
not look favorable when compared with the other potentials, 
the KS and the LDF one. 

We observe that the LDF potential is clearly the best choice 
among the cases studied. It will be employed for our further 
calculations in this work. The residual electron correlation 
can be estimated by taking the difference of the KS and LDF 
results. Table U shows that such difference yields a conserva- 
tive estimate in the case of 2p,j states, whereas for the 2s state 



it is only just consistent with the deviation of the LDF value 
from the exact result. We found that the choice of the param- 
eter x a = 1/2 in the definition of the KS potential [Eq. ( fTOl il 
yields a more conservative estimate for the 2s state. The resid- 
ual electron correlation was thus estimated as a difference of 
the LDF and KS results with the parameter x a = 1/2 for the 
2s state and x a = 2/3 for the 2pj states. The estimate ob- 
tained in this way decreases roughly as Z~ 2 with increase of 
Z and becomes entirely negligible (i.e., less than 5x 10 -6 a.u.) 
for Z > 21 for 2s state and Z > 30 for 2pj states. It can 
be mentioned that a small potential dependence of the MBPT 
results appears again for large values of Z and is due to an 
incomplete treatment of the negative-energy part of the Dirac 
spectrum. 

It is of interest to compare results of different calculations 
for lithium isoelectronic sequence in the low-Z region, where 
QED effects are small and do not significantly influence the 
comparison. In Figs. [2] and [3] we plot the absolute values of 
the difference of our MBPT values and the results of other 
calculations lT2ll l44l l45l l46l l47tl for the ionization energy of 
the 2s state and for the centroid of the 2p 1 / 2 and 2p^j 2 lev- 
els, respectively. Calculations by Chung 11451 l46l 14711 were 
carried out on nonrelativistic wave functions, first solving the 
Schrodinger equation and then adding the expectation value 
of the Breit-Pauli Hamiltonian. Relativistic calculations by 
Blundell et al. l44ll were performed within a combination of 
the coupled-cluster formalism and MBPT. They yielded rather 
accurate (by standards of MBPT) results but were carried out 
for lithium and beryllium only. The evaluation by Johnson et 
al. 1I2TI1 is a third-order MBPT calculation that stands most 
close to our treatment presented so far. In order to make 
a comparison with the results of other calculations, we sub- 
tracted the contribution of the mass-polarization operator H4 
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FIG. 3: The same as in Fig. [2] but for the center of gravity of the 
energies of the (ls) 2 2p 1 / 2 and (ls) 2 2p 3 / 2 states. The results by 
Chung are from 1 4711 . 

from the relativistic results of Refs. l45l l46l l47tl and the nu- 
clear recoil correction from those of Refs. 1 2jJ 44)] . 

We would like now to outline several differences between 
our present MBPT treatment and the one of Ref. [211] . First, 
we use a different starting point for the perturbative expansion 
(the LDF potential instead of the nonlocal DF potential). Sec- 
ond, we include the Breit part of the third-order energy cor- 
rection. This is an important issue since, as previously noted 
in Ref. [24], the third-order Coulomb contribution happens to 
be anomalously small, so that the corresponding Breit correc- 
tion is comparable with the Coulomb one even for low-Z sys- 
tems and becomes dominant in the high-Z region. Third, we 
include the negative-energy contribution for the corrections to 
the wave functions by summing over the complete Dirac spec- 
turm in Eqs. d2"TT > and (l22l . Fourth, we employ a much larger 
basis set than the one consisting of 20 positive-energy func- 
tions reported in Ref. 12111 . This becomes important for small 
values of Z, where significant numerical cancellations occur 
in our calculations of third-order energies. 

For the 2s state, we observe good agreement with Chung's 
values for the lowest- Z ions and a rapidly increasing devia- 
tion for higher- Z ions. The discrepancy is mainly due to his 
neglect of terms beyond the relative order (Za) 2 in the one- 
electron energies and in the one-photon exchange correction, 
which are included to all orders in relativistic calculations. In- 
deed, the leading correction beyond the Breit approximation 
can be written within the 1/Z expansion as 



SE{Z) = (Za) 6 a 60 + a(Za) 5 a 51 + 0[a 2 (Zaf 



(46) 



where the coefficient ago originates from the expansion of 
the Dirac energy, ago (2s) = —21/1024, and a^x comes 
from the one-photon exchange correction, 051 (2s) = 0.1175 
JH. Resulting contribution for, e.g., Z = 20 is SE(20) = 
—0.00266 a.u., which can be compared with the actual de- 
viation of our result from the Chung's one of —0.00275 a.u. 
The corresponding correction to the centroid of the 2p lev- 
els is much smaller (for Z = 10, it is about 1 x 10~ 6 a.u.), 



which explains much better agreement with Chung's results 
observed in Fig. [3] 

The comparison with other calculations for Z > 20 is post- 
poned to the next section where we supplement our MBPT 
values with a rigorous QED calculation of the two-photon ex- 
change correction. 



IV. QED TREATMENT 

In this section we present an ab initio QED treatment of the 
electron correlation up to second order of perturbation theory. 
Considering the first- and second-order corrections obtained 
in the previous section and given by Eqs. d43l and ( l44l . we ob- 
serve that the only part that involves some approximations is 
the two-photon exchange correction AE 1 -- 2 ' ^. Calculations 
of this correction for Li-like ions were previously performed 
in Refs. d [H [H [M [H • Our present task will be to carry 
out a similar calculation for the same model potentials as in 
our MBPT treatment. 

General formulas for the two-photon exchange correction 
are derived by the two-time Green's function method ll49Ll50fl . 
The correction can be conveniently represented by a sum of 
the two-electron ("2el") and three-electron ("3el") contribu- 
tions, each of which is subdivided into the irreducible ("ir") 
and the reducible ("red") parts, 



AE (2,0\2) = AE 2c\ + A£ ,2el + A£J : 



3el 



(47) 



The irreducible two-electron part is given by 



AE: 



2c 1 



E 

riiri 2 



I 

27 



i r iad,dir(w ) 7lin 2 ) 



(s c — u — ue ni ) (e v + u> - ue r . 

-Flad,ex(^,nin 2 ) 

(e v - uj - ue ni )(s c + w — ue n2 

^cr,dir(w,nin 2 ) 



(e c - lj - ue ni )(e v - uj - 
F CI<eK (uj, mn 2 ) 



1A£n2 ) 



(e v - u) - ue ni )(£v - w 



(48) 



where 



Flad.dir^, mn 2 ) = ^ Icvn in2 {u) In in2 cv(w) , (49) 

-Fiad,cx(w, nin 2 ) = - }^ I vcnin2 (uj) I nin2CV (u> - A) , 

(50) 

■Pcr,dir(<*>,rciTl2) = ^ Icn 2ni v(u) I nt vcn 2 M > ( 5l ) 



nin 2 l 



^2 En 2 n 1 v(^>) In 1 ccn 2 ( w — A) , 



(52) 



u= (1 — zO), A = e v — e c , and fix and /i 2 denote the 
momentum projections of the states n\ and n 2 , respectively. 
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The prime on the sum indicates that some intermediate states 
are excluded from the summation. Specifically, the omitted 
terms are: (e ni e„ 2 ) = {e c e v ), (e v e c ) in the first and second 
terms in the brackets, (e ni e„ 2 ) = (e c e v ) in the third term, and 
(e ni £„ 2 ) = (e c £ c ), (e v e v ) in the fourth term. The above con- 
ditions for omitting terms should be taken in the point-nucleus 
limit, i.e., the 2s and the 2f>]/ 2 states are treated in the same 
way despite the fact that their Dirac energies are separated by 
the finite nuclear size. 

The reducible two-electron part is 

/OO -i 

+2F CTyCX (-u>,vv) + 2F CIyCX (-uj + A vs , ss) 

--Flad,ex(w + A, cv) - -Flad,ox(-W + A, cv) 
--Flad,dir(w - A, VC) - Flad,dir(-W - A, vc) 
--Flad,ex(w, VC) - F lac L cx (-W, vc)} , (53) 

where s denotes the Dirac state whose energy is separated 
from the valence energy e v by the finite nuclear size only. 
(The corresponding contribution should be omitted if there is 
no such state, e.g., for v = 2p 3 / 2 .) 

The irreducible three-electron contribution reads 

AEt 1 = £M) P+Q E' 

PQ n 

Ip2 P3 n Q3 ( Ap3Q3 ) Ip\ „ Q\ Q2 ( Aqipi ) 

X I ' 

£ Ql + £ Q2 — £pi — £n 

(54) 

where the prime on the sum indicates that terms with the van- 
ishing denominator should be omitted. Finally, the reducible 
three-electron part is given by 

A-Efed = £ [^aa^(A) {I a b:ab ~ hv-bv) 

~\~2 ^avvb(^) Ibv-.av + ^ ^bvvai^) ^va-vb j 

(55) 

where we used the notation I ab . cd = I a bcd(A bd )-I bacd (A ad ), 
a and b denote the core electrons with the momentum projec- 
tions /i a and /!/,, n b = —^a, and v is a valence state with the 
momentum projection /j,y. 

We mention that the formulas d48l>-(f55T> reduce to the 
MBPT expression (f24t after neglecting (i) the energy depen- 
dence of the operator in the Coulomb gauge and (ii) the 
negative-energy part of the Dirac spectrum. Within this ap- 
proximation, all reducible parts vanish and the u> integration 
can be performed by Cauchy's theorem. 

The numerical evaluation of expressions d48b-(T55l> was per- 
formed by using the scheme described in detail in Refs. Ill ll 
[l2i[T5]l . Summations over the Dirac spectrum were carried out 
by using the dual-kinetic -balance basis set 1 39] constructed 
with B-splines; the typical value for the number of splines was 
65. The numerical uncertainty of our results was estimated to 
be of about 1 x 10~ 5 a.u. 



TABLE II: Electronic-structure contributions to the ionization po- 
tential of n = 2 states of Li-like Zn {Z — 30) for different model 
potentials, in a.u. 

Coul KS LDF 



2si/2 state: 



Dirac: 


-114.22811 


-102.26718 


-102.25381 


1 -photon: 


12.23022 


0.00568 


-0.00614 


2-photon(MBPT): 


-0.26803 


-0.00493 


-0.00644 


2-pnoton((jhl_>): 


—0.00012 


—0.00010 


—0.00010 


3-photon(MBPT): 


-0.00039 


0.00007 


0.00004 


Sum: 


— 102.26643 


— 102.26647 


— 102.26646 


2px/2 state: 








Dirac: 


-114.22864 


-99.85790 


-100.17793 


1 -photon: 


14.51407 


-0.26002 


0.05156 


2-photon(MBPT): 


-0.41102 


-0.00935 


-0.00066 


2-photon(QED): 


-0.00007 


-0.00006 


-0.00006 


3-photon(MBPT): 


-0.00141 


0.00016 


-0.00007 


Sum: 


-100.12706 


-100.12717 


-100.12717 


2p 3/2 state: 








Dirac: 


-112.83902 


-98.73159 


-99.05651 


1 -photon: 


14.18655 


-0.30055 


0.02088 


2-photon(MBPT): 


-0.38556 


-0.00819 


-0.00462 


2-photon(QED): 


0.00010 


0.00009 


0.00009 


3-photon(MBPT): 


-0.00209 


0.00011 


0.00003 


Sum: 


-99.04002 


-99.04013 


-99.04014 



The results of our calculations for the n = 2 states of Li- 
like zinc (Z = 30) and bismuth (Z = 83) are presented in 
Tables [TT] and [III] respectively. The entry "2-photon(QED)" 
represents the difference of the two-photon exchange correc- 
tion calculated within QED and within MBPT. For the mo- 
ment, we disregard uncertainties due to the nuclear size and 
the higher-order correlation; the purpose of Tables [Hi and [TTT1 
is to illustrate the convergence of results obtained for different 
model potentials. 

We observe that the total values obtained with the KS and 
the LDF potential almost coincide with each other for the both 
cases studied, whereas the Coulomb-potential value is slightly 
apart from the other two. The third-order correction turns out 
to be rather small for the KS and LDF potentials but is numer- 
ically significant in the Coulomb case. In the case of bismuth, 
we compare our results with the ones obtained previously by 
Sapirstein and Cheng [13], whose approach was very similar 
to ours. We report a good agreement with their calculation for 
the results obtained with the KS potential (a deviation for the 
2 s state is mainly due to a different nuclear radius used in that 
work). However, the results obtained with the Coulomb po- 
tential exhibit a small but noticeable difference. This is due to 
a certain disagreement in values for the three-photon exchange 
correction. For the KS potential and Z = 83, this correction 
is nearly negligible so that the discrepancy does not play any 
role. It can be mentioned that the differences between our 
Coulomb and KS results are much smaller than those reported 
by Sapirstein and Cheng. 

Our Coulomb results for the three-photon exchange correc- 
tion and Z = 30 presented in TableHTIare in reasonable agree- 
ment with Andreev et al. [14j], who obtained —0.00044 a.u. 
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TABLE III: Electronic-structure contributions to the ionization po- 
tential of n = 2 states of Li-like Bi (Z — 83) for different model 
potentials, in a.u. 





L.OU1 






■""1/2 SLdLC 








Dirac; 


— yo4.44io 


— y44.ZlOO 


— y44.0O0 / 


1 -photon; 


/in O/iQC 
4U.o4yo 


U.ZloU 


u.ooyz 






\J.U l-OO 


—0 01 




— U.UU1U 


n nm o 


n oooq 
— u.uuuy 


^ nhntnnfMRPTV 


n nm ^ 


n 0009 
u.uuuz 


n 0001 


oiim. 


— y44.Ulo / 


o/i /i m £7 
— y44.U10 / 


— y44.U!0o 


Ref. [131 


-944.0162 


-944.0155 




9n- >_ ctnt* 3 " 
£>\J\ / 2 ^ LtlLC. 








Dirac; 


ClQA S7SS 

— yo4.o / oo 


— yo4.o / 1 / 


— yoo.uouo 


1 -photon; 


CI QQ/1Q 

Ol.oo4o 


n fii on 


1 7COQ 


Z-pilOLOIl^lvlor 1 ). 


n 81 Oft 
— U.olUO 


n O/i 7*^ 


— u.uuoy 




n ooozt 

U.UUU4: 


000*} 


OOOzt 

U.UUUL: 


^-nhntnnnWRPTV 


U .\JVJOO 


000^ 


—0 0009 


ouni. 


QQ/1 OQQQ 


QO/1 OQQO 

— yo4.zyyz 


QO/1 OQQQ 

— yo4.zyyo 


Ref. ri3J 


-934.3019 


-934.2995 




2p 3/2 state: 








Dirac: 


-881.8298 


-839.8108 


-841.1109 


1 -photon: 


41.7248 


-0.7621 


0.5326 


2-photon(MBPT): 


-0.4958 


-0.0270 


-0.0216 


2-photon(QED): 


0.0106 


0.0102 


0.0102 


3-photon(MBPT): 


0.0008 


0.0000 


-0.0001 


Sum: 


-840.5893 


-840.5898 


-840.5898 


Ref. ri3J 


-840.5902 


-840.5896 





for the 2s state and —0.0015 a.u. for the 2p 1 / 2 state. How- 
ever, we disagree with them for the contribution induced by 
two Breit and one Coulomb interactions (which we calculated 
but did not include into the total values); our result for the 
2f>i/2 state is six times smaller. 

We performed our QED calculations of the two-photon ex- 
change correction for the nuclear charge numbers in the range 
10 < Z < 92. At the lowest value considered, Z = 10, the 
difference of the QED and MBPT results is already smaller 
than the accuracy we are presently interested in, 1 x 10~ 5 a.u., 
for all n = 2 states. The difference scales as Z 3 and is con- 
sidered negligible for Z < 10. In Fig.|4j we plot this differ- 
ence scaled by a factor of (Za) 3 as a function of the nuclear 
charge number. We observe that the "pure" QED part of the 
two-photon exchange correction is remarkably small for the 
2s and 2pi/ 2 states even in the high-Z region but is much 
larger in the case of the 2p 3 / 2 state. 

Our total results for the electronic-structure corrections to 
the ionization potential of the n — 2 states of Li-like ions 
are collected in Table [VTJ The values for the lightest atoms 
with Z = 3 and 4 are given only for the completeness; by 
far more accurate calculations are available for these systems 
||42l|43i|44ll. 

The column labeled "Dirac" contains the energy values (mi- 
nus mc 2 ) obtained from the Dirac equation with the Fermi- 
like nuclear potential. The values for the the nuclear-charge 
root-mean-square (rms) radii and their uncertainties are listed 
in the second column of the table. They were taken from 



Ref. |[5lll except for few cases with no experimental data 
available (Z = 43, 61, 85, 89, and 91), in which case we 
used the interpolation formula from Ref. ll52tl and assigned 
an uncertainty of 1% to these values. The dependence of 
the Dirac value on the nuclear model was conservatively esti- 
mated by comparing the results obtained within the Fermi and 
the homogeneously-charged-sphere models. 

The entries "1 -photon", "2-ph.[MBPT]", and "3-photon" 
contain results for the one-, two-, and three-photon exchange 
corrections evaluated within MBPT with the LDF potential 
and given by Eqs. d43l)-(l45b. respectively. The entry "2- 
ph.[QED]" represents the difference of the two-photon ex- 
change correction evaluated within QED and MBPT, i.e., the 
difference of Eqs. (f4Tb and d44T >. The numerical uncertainty 
of 1 x 10~ 5 a.u. is not specified explicitly in Table PVll but 
included into the total error estimate. 

The column "h.o." contains errors due the higher-order ef- 
fects neglected in the present investigation. Our estimations 
for these effects consist of two parts that are added quadrat- 
ically, the residual electron correlation and the QED part of 
the three-photon exchange correction. The residual correla- 
tion was discussed in the previous section; it has its maximum 
for Z = 3 and decreases rapidly when Z increases. Mak- 
ing a comparison with more precise calculations available for 
lithium and beryllium allow us to be reasonably confident of 
this part of our error estimate. 

An estimation of residual three-photon QED effects is more 
difficult to make reliably. Fig.|4]leads us to surmise that QED 
part of the two-photon exchange correction is anomalously 
small for the 2s and 2pi/ 2 states; moreover, it changes its sign 
when Z is varied (for 2p states). An estimate based on the 
ratio of the QED and MBPT two-photon contributions would 
thus likely to underestimate the three-photon QED effects. We 
choose to base our estimation on the three-photon MBPT con- 
tribution induced by two Breit and one Coulomb interactions; 
let us denote it by AEI A q PT (B x B). It has the same scal- 
ing order as the QED contribution [a 3 (Za) 2 ] and it does not 
change its sign through the Z range of interest. For the 2s and 
2pi/ 2 states, we take the absolute value of this contribution 
for the estimation of three-photon QED effects. For the 2p 3 / 2 
state, we multiply its value by the ratio 



A-EmbptO 8 x b ) 

whose numerical contribution varies from 1 in the middle- Z 
region to 5 in the high-Z region. 

We are now in a position to compare our numerical results 
for the electronic-structure part of the 2pj-2s transition ener- 
gies with results of other investigations. We selected three ex- 
tensive investigations of the electronic structure of Li-like ions 
to compare with, which were accomplished by three indepen- 
dent methods: MBPT 12111 . the multiconfigurational Dirac - 
Fock (MCDF) method 02211 . and the configuration interaction 
(CI) method 12311 . It should be noted that none of these cal- 
culations accounted for the QED part of the two-photon ex- 
change correction, which is included into consideration in the 
present work. 
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FIG. 4: The difference of the QED and MBPT results for the two-photon exchange correction for the 2s, 2p\/2, and 2p 3 / 2 states of Li-like 
ions, scaled by a factor of (Za) 3 . 
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FIG. 5: The deviation of the present results for the the electronic-structure part of the 2pi/ 2 -2a and 2p 3 / 2 -2s transition energies from the 
previous calculations (dots, MBPT [21.]; stars, CI (23ll ; triangles, MCDF 1 22]). The dotted line represents the estimated error of our results due 
to the residual electron correlation; the dashed line denotes the QED part of the two-photon exchange correction accounted for in the present 
investigation. 
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In Fig. [3] we plot the deviation of our numerical values 
for the electronic-structure part of the 2p 3 / 2 -2s and 2p 1 / 2 -2s 
transition energies of Li-like ions from the results of previ- 
ous calculations. The electronic-structure part can be unam- 
biguously isolated from the data presented in Refs. I21I l23ll ; 
from Ref. fl22ll . we had to subtract the mass-polarization term, 
which is explicitly given there for Z — 15, 54 and 92 only. In 
Fig. 13 the dotted line indicates the estimated error of our re- 
sults due to the residual correlation; the dashed line stands for 
the QED part of the two-photon exchange correction, which 
is accounted for in our calculation but is omitted in previous 
studies. 

We observe a distinct and nearly Z-independent deviation 
of our results from the MCDF values, which is on the level 
of about 0.005-0.01 a.u. for the 2pi/ 2 -2s transition and of 
about 0.003 a.u. for the 2p 3 / 2 -2s transition. Agreement with 
the MBPT and CI results is much better; one can observe that 
the main part of the deviation is, as expected, due to the QED 
effects not accounted for in the previous studies. 

V. TRANSITION ENERGIES IN LI-LIKE IONS 

A number of important corrections should be added to the 
electronic-structure part of transition energies addressed to in 
the previous section in order to allow an adequate compari- 
son with experimental data. We now briefly discuss various 
theoretical contributions to the 2p 1 / 2 -2s and 2p 3 / 2 -2s transi- 
tion energies in Li-like ions. Our discussion is summarized 
by Tables [IV] and [V] where various individual contributions 
are collected for several medium-Z ions. 

The largest effect to be added to the electronic-structure 
part is the one-loop QED contribution. It consists of the self- 
energy and vacuum-polarization corrections and is presently 
well-understood; we refer the reader to the review lfl8ll for 
the details. The numerical values listed in Tables ITVl and M 
represent the one-loop QED correction calculated on the hy- 
drogenic wave functions. 

The next effect is the screening of the one-loop QED cor- 
rections by other electrons. Rigorous evaluations of the first- 
order (in 1/Z) screening effects were performed for Li-like 
10ns in our previous investigations 

HEM; 

a similar cal- 
culation was carried out for Li-like bismuth by Sapirstein and 
Cheng [13]. Higher-order screening effects have not been cal- 
culated up to now. The error due to their neglect has to be 
estimated with some care since it presently yields the dom- 
inant theoretical uncertainty for medium-Z ions. To obtain 
such an estimate, we recall that, within the Za expansion, the 
dominant part of the screening of one-electron QED correc- 
tions can be described ll58l 15911 by incorporating the correct 
electron density at the nucleus into the hydrogenic formulas. 
For the 2p-2s transition, the nonrelativistic electron density at 
the origin is given by [60] 

(E< <K r *))(i s )2 2p - (Ei S (^))(l S )2 2s 

_ GM_ r_i/8 + 0.339 893 Z- 1 - 0.139 Z~ 2 + 0(Z- 3 )] . 

TT 

(56) 



The ratio of the second and the first term in the above ex- 
pression (1=3 — 2.7 /Z) is very close to the actual ratio of the 
first-order screening contribution to the hydrogenic one. We 
estimate the ratio of the higher-order effects to the first-order 
screening by taking the ratio of the Z~ 2 term in Eq. (|56| > to 
the Z -1 one and multiplying it by a conservative factor of 2, 
the resulting scaling factor thus being 0.8/ Z. 

Another important effect to be taken into account is the 
recoil correction. Rigorous QED calculations of the leading 
term of the 1/Z expansion of this correction were performed 
in Refs. lfl9l I20I1 . The higher-order (in 1/Z) part of the recoil 
effect has been addressed to in the case of Li-like ions only 
nonrelativistic ally up to now. Within the nonrelativistic ap- 
proximation, we obtain it by evaluating the reduced-mass cor- 
rection to the electronic-structure part of order 1/Z and higher 
and (for the 2p states) by adding the part of the mass polariza- 
tion of order 1/Z and higher inferred from the Hughes-Eckart 
formula l6lll . A 100% uncertainty was acribed to the part of 
the recoil effect obtained within the nonrelativistic approxi- 
mation. 

Finally, we should account for the two-loop QED effects. 
The first complete all-order calculation of the two-loop one- 
electron QED corrections for n = 2 states was recently ac- 
complished in Ref. lfl7ll for several ions with Z > 60. In order 
to complete our compilation for ions with Z = 26, 28, 36, and 
47 in Tables ITVl and IVl we need values for the two-loop QED 
corrections outside the Z range covered in Ref. B17I1 . To this 
end, we performed calculations of the two-loop diagrams in- 
volving closed electron loops. In notations of Ref. [17], these 
are the subsets "SEVP", "VPVP", and "S(VP)E". The results 
are presented in Tables [TVl and [VI the error bars indicated are 
due to the free-loop approximation employed in the evaluation 
of the VPVP and S(VP)E subsets, see Ref. [fTvb for details. 

The remaining two-loop self-energy correction (the 
"SESE" subset) is much more difficult to calculate and we ob- 
tain its numerical values by an extrapolation. For the 2s state, 
the extrapolation was performed in two steps. First, we obtain 
numerical values for the Is state, which was done by interpo- 
lating the numerical results of Ref. 116211 . Second, we obtain 
results for the weighted difference of the 2s and Is correc- 
tions, A s = 8 SE 2s — 8Ei s . This was achieved by subtracting 
all Za-expansion contributions known (see Refs. B63l 64 11 and 
references therein) from the all-order results of Ref. jl7ll and 
extrapolating the higher-order remainder towards the values 
of Z of interest. An uncertainty of 30% was ascribed to these 
results. For the 2p states, the correction is much smaller and, 
for our purposes, it is sufficient to obtain just the boundaries 
for the higher-order remainder. The comparison [65}] of the 
all-order numerical values with the Za-expansion results sug- 
gests an estimate of ±2 a 2 (Za) 6 /(8 ir 2 ) for the higher-order 
remainder. 

Considering the data presented in Tables [TVl and [VI we no- 
tice a remarkably good agreement of our total values with 
the experimental results. The theoretical accuracy is lower 
than the experimental one in the case of iron and nickel but 
significantly better for krypton and silver. It should be men- 
tioned that the leading theoretical uncertainty for medium-Z 
ions stems presently from the higher-order screening and re- 
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TABLE IV: Individual contributions to the 2p!/ 2 -2s transition energy in Li-like iron, nickel, and krypton, in eV (1 a.u. = 27.211 383 eV). 





Subset 


Z = 26 


Z = 28 


Z = 36 


Electronic structure 




49.1031 (5) 


53.6036(3) 


72.8013 (6) 


One-loop QED 




-0.5565 


-0.7168 


-1.6859 


Screened QED 




0.0649 (20) 


0.0778 (22) 


0.1433 (32) 


Recoil 




-0.0121(19) 


-0.0137(20) 


-0.0163(18) 


Two-loop QED 


SESE 


0.0004(1) 


0.0006(1) 


0.0022 (4) 




SEVP 


-0.0002 


-0.0002 


-0.0008 




VPVP 


0.0004 


0.0005 


0.0015(1) 




S(VP)E 


-0.0001 


-0.0001 


-0.0002 (2) 


Total theory 




48.6000 (28) 


52.9517(30) 


71.2451 (37) 


Experiment 




48.5997(10) [53] 


52.9501 (11) [54] 


71.243 (8) [4] 






48.6033(19) [55J 


52.9497(22) [55J 


71.241 (11) [55J 






48.6001 (19) [56J 


52.947 (4) [57] 





TABLE V: Individual contributions to the 2p3/2-2s transition energy in Li-like iron, nickel, and silver, in eV. 





Subset 


Z = 26 


Z = 28 


Z = 47 


Electronic structure 




65.0333 (5) 


75.5634 (3) 


307.1988(16) 


One-loop QED 




-0.5119 


-0.6574 


-3.7439 


Screened QED 




0.0551 (17) 


0.0658 (19) 


0.2316 (39) 


Recoil 




-0.0123 (20) 


-0.0140 (21) 


-0.0243(19) 


Two-loop QED 


SESE 


0.0003 (1) 


0.0005 (1) 


0.0072 (20) 




SEVP 


-0.0002 


-0.0002 


-0.0032 




VPVP 


0.0004 


0.0005 


0.0047 (5) 




S(VP)E 


-0.0001 


-0.0001 


-0.0007(7) 


Total theory 




64.5647 (26) 


74.9585 (28) 


303.6704 (52) 


Experiment 




64.5656 (17) [53] 


74.9602 (22) [54J 


303.67 (3) [2J 






64.571 (7) [55J 


74.962 (4) [55J 








64.560 (3) \56} 


74.958 (7) [57] 





coil effects. In this respect, the situation for medium-Z ions 
is different from the one encountered in the high-Z region, 
where the uncertainties due to the finite nuclear size effect 
and due to the two-loop QED corrections become prominent 

The uncertainty due to the screening of QED effects can be 
reduced by calculating the one-loop and the first-order screen- 
ing QED corrections not on hydrogenic wave functions but in 
the presence of a local screening potential, as it was done in 
the case of bismuth in Ref . lfl3ll . The uncertainty of the recoil 
effect can also be improved by calculating the higher-order 
(in 1/Z) recoil correction within the leading relativistic ap- 
proximation. This means that the theoretical accuracy can be 
pushed even further in the near future. 



For the first two terms of the expansion, rigorous QED calcu- 
lations were performed, whereas the third-order contribution 
was evaluated within the MBPT approximation. Errors due to 
truncation of the perturbative expansion and due to usage of 
the MBPT approximation in the third-order correction were 
estimated. 

Collecting all theoretical contributions available for the 
2p,/-2s transition energies, we observed good agreement with 
experimental results for medium-Z ions. The dominant un- 
certainty of the theoretical results in the medium-Z region is 
shown to originate from the higher-order screening of QED 
corrections and from the recoil effect. Further improvement 
of the theoretical accuracy is anticipated for medium-Z ions. 



VI. CONCLUSION 

In this paper we have presented a systematic QED treatment 
of the electron correlation for n = 2 states of Li-like ions. The 
treatment relies on the perturbative expansion, with a local 
model potential included into the zeroth-order approximation. 
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TABLE VI: Electronic-structure contributions to the ionization potential of the 



2 states of Li-like ions, in a.u. 



~Z <r I > m 



State 



Dirac 


1 -photon 


2ph.[MBPT] 2ph.[QED] 


3-photon 




h.o. 


Total 




-0 


19601 


-0 


00031 


-0 


00184 





00000 





00001 


± 





00038 


-0 


19815 


(38) 


-0 


12867 


-0 


00109 


-0 


00045 





00000 


-0 


00009 


± 





00039 


-0 


13029 


(39) 


-0 


12867 


-0 


00108 


-0 


00046 





00000 


-0 


00009 


± 





00039 


-0 


13029 


(39) 


-0 


66520 


-0 


00117 


-0 


00293 





00000 





00000 


± 





00021 


-0 


66930 


(21) 


-0 


51971 


-0 


00303 


-0 


00095 





00000 


-0 


00016 


± 





00018 


-0 


52386 


(18) 


-0 


51967 


-0 


00298 


-0 


00102 





00000 


-0 


00016 


± 





00018 


-0 


52383 


(18) 


-1 


38848 


-0 


00202 


-0 


00352 





00000 





00000 


± 





00012 


-1 


39402 


(12) 


-1 


16794 


-0 


00424 


-0 


00133 





00000 


-0 


00015 


± 





00011 


-1 


17366 


(12) 


-1 


16775 


-0 


00416 


-0 


00147 





00000 


-0 


00013 


± 





00011 


-1 


17351 


(12) 


-2 


36360 


-0 


00311 


-0 


00354 





00000 


-0 


00003 


± 





00008 


-2 


37027 


(8) 


-2 


06991 


-0 


00472 


-0 


00160 





00000 


-0 


00012 


± 





00006 


-2 


07636 


(6) 


-2 


06933 


-0 


00462 


-0 


00182 





00000 


-0 


00010 


± 





00006 


-2 


07588 


(6) 


-3 


58947 


-0 


00416 


-0 


00400 





00000 


-0 


00001 


± 





00005 


-3 


59764 


(5) 


-3 


22370 


-0 


00526 


-0 


00176 





00000 


-0 


00011 


± 





00006 


-3 


23083 


(6) 


-3 


22234 


-0 


00516 


-0 


00208 





00000 


-0 


00008 


± 





00006 


-3 


22966 


(6) 


-5 


06676 


-0 


00508 


-0 


00428 





00000 





00000 


± 





00003 


-5 


07613 


(3) 


-4 


62925 


-0 


00550 


-0 


00186 





00000 


-0 


00010 


± 





00005 


-4 


63671 


(5) 


-4 


62651 


-0 


00544 


-0 


00229 





00000 


-0 


00006 


± 





00005 


-4 


63430 


(5) 


-6 


79562 


-0 


00599 


-0 


00432 





00000 


-0 


00001 


± 





00003 


-6 


80593 


(3) 


-6 


28658 


-0 


00539 


-0 


00192 





00000 


-0 


00008 


± 





00003 


-6 


29398 


(4) 


-6 


28163 


-0 


00541 


-0 


00247 





00000 


-0 


00005 


± 





00003 


-6 


28955 


(4) 


-8 


77583 


-0 


00696 


-0 


00457 





00000 





00000 


± 





00003 


-8 


78736 


(3) 


-8 


19547 


-0 


00531 


-0 


00194 





00000 


-0 


00008 


± 





00003 


-8 


20281 


(3) 


-8 


18719 


-0 


00547 


-0 


00261 





00000 


-0 


00003 


± 





00003 


-8 


19531 


(3) 


-11 


00800 


-0 


00813 


-0 


00465 





00000 





00000 


± 





00002 


-11 


02078 


(3) 


-10 


35648 


-0 


00496 


-0 


00194 


-0 


00001 


-0 


00007 


± 





00002 


-10 


36347 


(3) 


-10 


34345 


-0 


00532 


-0 


00274 





00000 


-0 


00003 


± 





00002 


-10 


35154 


(3) 


-13 


49330 


-0 


00863 


-0 


00470 





00000 





00000 


± 





00002 


-13 


50663 


(2) 


-12 


76996 


-0 


00437 


-0 


00192 


-0 


00001 


-0 


00007 


± 





00002 


-12 


77633 


(2) 


-12 


75036 


-0 


00500 


-0 


00286 





00000 


-0 


00002 


± 





00002 


-12 


75825 


(2) 


-16 


23067 


-0 


00999 


-0 


00476 


-0 


00001 





00001 


± 





00002 


-16 


24542 


(2) 


-15 


43621 


-0 


00367 


-0 


00189 


-0 


00001 


-0 


00007 


± 





00002 


-15 


44184 


(2) 


-15 


40784 


-0 


00466 


-0 


00297 





00000 


-0 


00001 


± 





00002 


-15 


41549 


(2) 


-19 


22244 


-0 


01040 


-0 


00488 


-0 


00001 





00001 


± 





00001 


-19 


23772 


(2) 


-18 


35582 


-0 


00277 


-0 


00184 


-0 


00001 


-0 


00007 


± 





00002 


-18 


36050 


(2) 


-18 


31603 


-0 


00422 


-0 


00307 





00000 


-0 


00001 


± 





00002 


-18 


32333 


(2) 


-22 


46840 


-0 


01073 


-0 


00500 


-0 


00001 





00001 


± 





00001 


-22 


48413 


(2) 


-21 


52936 


-0 


00166 


-0 


00178 


-0 


00001 


-0 


00006 


± 





00001 


-21 


53287 


(2) 


-21 


47501 


-0 


00367 


-0 


00317 





00000 





00000 


± 





00001 


-21 


48185 


(2) 


-25 


96925 


-0 


01096 


-0 


00511 


-0 


00001 





00001 


± 





00001 


-25 


98532 


(2) 


-24 


95748 


-0 


00031 


-0 


00172 


-0 


00001 


-0 


00006 


± 





00001 


-24 


95959 


(!) 


-24 


88489 


-0 


00301 


-0 


00326 





00000 





00000 


± 





00001 


-24 


89115 


(!) 


-29 


72572 


-0 


01110 


-0 


00521 


-0 


00002 





00001 


± 





00001 


-29 


74203 


(!) 


-28 


64087 





00126 


-0 


00165 


-0 


00001 


-0 


00006 


± 





00001 


-28 


64133 


(!) 


-28 


54578 


-0 


00224 


-0 


00335 





00000 





00001 


± 





00001 


-28 


55136 


(!) 


-33 


73717 


-0 


01263 


-0 


00523 


-0 


00002 





00001 


± 





00001 


-33 


75504 


(!) 


-32 


58029 





00310 


-0 


00157 


-0 


00001 


-0 


00006 


± 





00001 


-32 


57884 


(!) 


-32 


45782 


-0 


00135 


-0 


00344 





00000 





00001 


± 





00001 


-32 


46260 


(!) 


-38 


00721 


-0 


01264 


-0 


00533 


-0 


00002 





00002 


± 





00001 


-38 


02519 


(!) 


-36 


77654 





00521 


-0 


00150 


-0 


00002 


-0 


00006 


± 





00001 


-36 


77291 


(!) 


-36 


62116 


-0 


00033 


-0 


00353 





00000 





00001 


± 





00001 


-36 


62500 


(!) 


-42 


53541 


-0 


01252 


-0 


00544 


-0 


00002 


o 


00002 


± 


o 


00001 


-42 


55337 


(!) 


-41 


23051 





00762 


-0 


00142 


-0 


00002 


-0 


00006 


± 





00001 


-41 


22439 


(!) 


-41 


03595 





00082 


-0 


00362 





00001 





00001 


± 





00001 


-41 


03873 


(!) 


-47 


32269 


-0 


01227 


-0 


00557 


-0 


00003 





00002 


± 





00001 


-47 


34053 


(!) 


-45 


94308 





01031 


-0 


00134 


-0 


00002 


-0 


00006 


± 





00001 


-45 


93420 


(!) 


-45 


70234 





00209 


-0 


00371 





00001 





00002 


± 





00001 


-45 


70394 


(!) 


-52 


37013 


-0 


01189 


-0 


00567 


-0 


00003 





00002 


± 





00001 


-52 


38770 


(!) 


-50 


91531 





01335 


-0 


00126 


-0 


00003 


-0 


00006 


± 





00001 


-50 


90331 


(!) 


-50 


62056 





00352 


-0 


00380 





00001 





00002 


± 





00001 


-50 


62081 


(1) 



3 2.431 (28) 

4 2.518(11) 

5 2.406 (29) 

6 2.470 (2) 

7 2.558(7) 

8 2.701 (6) 

9 2.898(2) 

10 3.005(2) 

11 2.994(2) 

12 3.057(2) 

13 3.061 (4) 

14 3.122(2) 

15 3.189(2) 

16 3.261 (2) 

17 3.365(15) 

18 3.427(2) 

19 3.435 (2) 

20 3.476 (1) 

21 3.544 (2) 

22 3.591 (2) 



2S1/2 
2P1/2 

2p 3 /2 

2S1/2 
2P1/2 
2P3/2 

2si/ 2 

2P1/2 
2P3/2 
2*1/2 
2P1/2 
2P3/2 
2Si/ 2 
2P1/2 
2P3/2 
2*1/2 

2pi/ 2 

2P3/2 
2Si/ 2 
2P1/2 
2P3/2 
2Sl/ 2 
2P1/2 
2P3/2 
2Si/ 2 
2P1/2 
2P3/2 

2si/ 2 

2P1/2 
2^3/2 

2si/ 2 

2P1/2 
2^3/2 
2Si/ 2 

2pi/ 2 

2^3/2 

2si/ 2 

2P1/2 
2P3/2 
2Si/ 2 
2P1/2 
2P3/2 
2Si/ 2 
2f>l/2 
2^3/2 
2Si/ 2 
2J51/2 
2P3/2 
2Si/ 2 
2P1/2 
2P3/2 
2Sl/ 2 
2j5l/2 
2f>3/2 
2Si/ 2 

2pi/ 2 

2P3/2 
2Si/ 2 

2pi/ 2 

2^3/2 
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TABLE VI: Electronic-structure contributions (continued) 



z 




State 


Dirac 


1 -photon 


2ph.[MBPT] 2ph.[QED] 


3-photon 


h.o. 


Total 




23 


3.599 (2) 


2Si/ 2 


-57 


67878 


-0.01138 


-0.00577 


-0.00004 


0.00003 


± 0.00001 


-57 


69595 


(!) 






2P1/2 


-56 


14820 


0.01673 


-0.00118 


-0.00003 


-0.00007 


± 0.00001 


-56 


13275 


(!) 






2P3/2 


-55 


79077 


0.00510 


-0.00389 


0.00002 


0.00002 


± 0.00001 


-55 


78953 


(!) 


24 


3.642 (2) 


2*1/2 


-63 


24979 


-0.01070 


-0.00587 


-0.00005 


0.00003 


± 0.00001 


-63 


26638 


(!) 






2pi/ 2 


-61 


64287 


0.02048 


-0.00110 


-0.00004 


-0.00007 


± 0.00001 


-61 


62359 


(!) 






2p 3 /2 


-61 


21319 


0.00684 


-0.00399 


0.00002 


0.00002 


± 0.00001 


-61 


21029 


(!) 


25 


3.706 (2) 


2si/ 2 


-69 


08434 


-0.00989 


-0.00598 


-0.00006 


0.00003 


± 0.00001 


-69 


10023 


(!) 




2pi/ 2 


-67 


40049 


0.02459 


-0.00102 


-0.00004 


-0.00007 


± 0.00001 


-67 


37702 


(!) 






2^3/2 


-66 


88801 


0.00874 


-0.00409 


0.00003 


0.00002 


± 0.00001 


-66 


88330 


(!) 


26 


3.737(2) 


2si/ 2 


-75 


18369 


-0.00890 


-0.00609 


-0.00006 


0.00003 


± 0.00001 


-75 


19871 


(!) 




2pi/2 


-73 


42227 


0.02912 


-0.00094 


-0.00004 


-0.00007 


± 0.00001 


-73 


39421 


(!) 






2p 3 /2 


-72 


81547 


0.01081 


-0.00419 


0.00004 


0.00002 


± 0.00001 


-72 


80878 


(!) 


27 


3.788 (2) 


2si/ 2 


-81 


54918 


-0.00775 


-0.00620 


-0.00007 


0.00003 


± 0.00001 


-81 


56317 


(!) 






2Pl/2 


-79 


70952 


0.03405 


-0.00087 


-0.00005 


-0.00007 


± 0.00001 


-79 


67645 


(!) 






2p 3 /2 


-78 


99579 


0.01305 


-0.00429 


0.00005 


0.00002 


± 0.00001 


-78 


98696 


(!) 


28 


3.775(1) 


2Si/ 2 


-88 


17955 


-0.00917 


-0.00621 


-0.00008 


0.00003 


± 0.00001 


-88 


19497 


(!) 






2pi/ 2 


-86 


26358 


0.03942 


-0.00080 


-0.00005 


-0.00007 


± 0.00001 


-86 


22508 


(!) 






2P3/2 


-85 


42924 


0.01548 


-0.00440 


0.00006 


0.00003 


± 0.00001 


-85 


41807 


(!) 


29 


3.882 (2) 


2si/ 2 


-95 


08143 


-0.00775 


-0.00632 


-0.00009 


0.00003 


± 0.00001 


-95 


09556 


(!) 






2pi/ 2 


-93 


08590 


0.04525 


-0.00073 


-0.00005 


-0.00007 


± 0.00001 


-93 


04150 


(!) 






2P3/2 


-92 


11605 


0.01808 


-0.00451 


0.00007 


0.00003 


± 0.00001 


-92 


10238 


(!) 


30 


3.929(1) 


2Si/ 2 


-102 


25381 


-0.00614 


-0.00644 


-0.00010 


0.00004 


± 0.00001 


-102 


26646 


(!) 




2pi/ 2 


-100 


17793 


0.05156 


-0.00066 


-0.00006 


-0.00007 


± 0.00001 


-100 


12717 


(!) 






2p 3 /2 


-99 


05651 


0.02088 


-0.00462 


0.00009 


0.00003 


± 0.00001 


-99 


04014 


(!) 


31 


3.997(2) 


2Si/ 2 


-109 


69828 


-0.00434 


-0.00656 


-0.00011 


0.00004 


± 0.00001 


-109 


70926 


(!) 






2pi/ 2 


-107 


54126 


0.05835 


-0.00060 


-0.00006 


-0.00007 


± 0.00001 


-107 


48365 


(!) 






2P3/2 


-106 


25089 


0.02388 


-0.00474 


0.00010 


0.00003 


± 0.00001 


-106 


23162 


(2) 


32 


4.074(1) 


2si/ 2 


-117 


41648 


-0.00234 


-0.00668 


-0.00013 


0.00004 


± 0.00001 


-117 


42559 


(!) 






2Pl/2 


-115 


17749 


0.06565 


-0.00055 


-0.00007 


-0.00007 


± 0.00001 


-115 


11252 


(2) 






2f>3/2 


-113 


69948 


0.02708 


-0.00486 


0.00012 


0.00003 


± 0.00001 


-113 


67710 


(2) 


33 


4.097(2) 


2Si/ 2 


-125 


41020 


-0.00009 


-0.00681 


-0.00014 


0.00004 


± 0.00001 


-125 


41720 


(!) 






2P1/2 


-123 


08834 


0.07352 


-0.00049 


-0.00007 


-0.00007 


± 0.00001 


-123 


01546 


(2) 






2^3/2 


-121 


40257 


0.03049 


-0.00499 


0.00014 


0.00003 


± 0.00001 


-121 


37690 


(2) 


34 


4.140 (2) 


2si/ 2 


-133 


68119 


0.00236 


-0.00694 


-0.00015 


0.00004 


± 0.00001 


-133 


68588 


(2) 




2pi/ 2 


-131 


27557 


0.08193 


-0.00044 


-0.00007 


-0.00008 


± 0.00001 


-131 


19424 


(2) 






2P3/2 


-129 


36048 


0.03411 


-0.00512 


0.00017 


0.00003 


± 0.00002 


-129 


33129 


(2) 


35 


4.163 (2) 


2si/ 2 


-142 


23136 


0.00503 


-0.00707 


-0.00017 


0.00004 


± 0.00001 


-142 


23352 


(2) 






2P1/2 


-139 


74102 


0.09090 


-0.00040 


-0.00008 


-0.00008 


± 0.00002 


-139 


65067 


(2) 






2^3/2 


-137 


57354 


0.03795 


-0.00525 


0.00019 


0.00003 


± 0.00002 


-137 


54062 


(2) 


36 


4.188(1) 


2Si/ 2 


-151 


06266 


0.00794 


-0.00720 


-0.00018 


0.00004 


± 0.00001 


-151 


06206 


(2) 






2P1/2 


-148 


48661 


0.10047 


-0.00036 


-0.00008 


-0.00008 


± 0.00002 


-148 


38666 


(2) 






2P3/2 


-146 


04206 


0.04201 


-0.00539 


0.00022 


0.00003 


± 0.00002 


-146 


00519 


(2) 


37 


4.203 (2) 


2si/ 2 


-160 


17714 


0.01110 


-0.00734 


-0.00020 


0.00005 


± 0.00001 


-160 


17353 


(2) 






2pi/ 2 


-157 


51436 


0.11067 


-0.00032 


-0.00009 


-0.00008 


± 0.00002 


-157 


40418 


(2) 






2f>3/2 


-154 


76641 


0.04630 


-0.00554 


0.00026 


0.00003 


± 0.00002 


-154 


72536 


(2) 


38 


4.220(1) 


2Si/ 2 


-169 


57690 


0.01450 


-0.00748 


-0.00021 


0.00005 


± 0.00001 


-169 


57004 


(2) 




2^1/2 


-166 


82632 


0.12149 


-0.00029 


-0.00009 


-0.00008 


± 0.00002 


-166 


70529 


(2) 






2P3/2 


-163 


74693 


0.05083 


-0.00569 


0.00029 


0.00003 


± 0.00002 


-163 


70147 


(3) 


39 


4.242 (2) 


2si/ 2 


-179 


26414 


0.01817 


-0.00762 


-0.00023 


0.00005 


± 0.00002 


-179 


25377 


(2) 




2P1/2 


-176 


42466 


0.13299 


-0.00026 


-0.00010 


-0.00008 


± 0.00002 


-176 


29211 


(2) 






2P3/2 


-172 


98399 


0.05559 


-0.00584 


0.00033 


0.00003 


± 0.00003 


-172 


93388 


(3) 


40 


4.270(1) 


2si/ 2 


-189 


24117 


0.02212 


-0.00777 


-0.00025 


0.00005 


± 0.00002 


-189 


22702 


(2) 






2pi/ 2 


-186 


31165 


0.14518 


-0.00024 


-0.00010 


-0.00009 


± 0.00002 


-186 


16690 


(2) 






4P3/2 


-182 


47797 


U.UUUUU 


-0.00600 


0.00037 


n nnnn* 1 ! 

u.uuuuo 


zn u.uuuuo 


-182 


42297 


O j 


41 


4.324(2) 


2Si/ 2 


-199 


51035 


0.02636 


-0.00792 


-0.00026 


0.00005 


± 0.00002 


-199 


49212 


(2) 






2P1/2 


-196 


48962 


0.15808 


-0.00023 


-0.00010 


-0.00009 


± 0.00002 


-196 


33196 


(3) 






2P3/2 


-192 


22925 


0.06585 


-0.00617 


0.00042 


0.00003 


± 0.00003 


-192 


16913 


(3) 


42 


4.409(1) 


2si/ 2 


-210 


07409 


0.03089 


-0.00808 


-0.00028 


0.00005 


± 0.00002 


-210 


05151 


(2) 




2pi/ 2 


-206 


96097 


0.17172 


-0.00022 


-0.00011 


-0.00009 


± 0.00003 


-206 


78967 


(3) 






2P3/2 


-202 


23825 


0.07135 


-0.00634 


0.00047 


0.00003 


± 0.00003 


-202 


17274 


(4) 



TABLE VI: Electronic-structure contributions (continued) 



z 


■} 1 l*> 

< r v /2 


State 


Dirac 


1 -photon 


2ph.[MBPT] 2ph.[QED] 


3-photon 


h.o. 


Total 




43 


4.424 (44) 


2Si/ 2 


-220 


93510 (8) 





03569 


-0.00823 


-0.00030 


0.00005 


± 0.00002 


-220 


90789 


(8) 






2P1/2 


-217 


72825 





18610 


-0.00022 


-0.00011 


-0.00009 


± 0.00003 


-217 


54257 


(3) 






2p 3 /2 


-212 


50536 





07711 


-0.00652 


0.00053 


0.00003 


± 0.00004 


-212 


43422 


(4) 


44 


4.482 (2) 


2S1/2 


-232 


09593 





04086 


-0.00840 


-0.00032 


0.00006 


± 0.00002 


-232 


06372 


(2) 






2P1/2 


-228 


79408 





20132 


-0.00022 


-0.00011 


-0.00009 


± 0.00003 


-228 


59319 


(3) 






2P3/2 


-223 


03102 





08313 


-0.00671 


0.00059 


0.00002 


± 0.00004 


-222 


95399 


(4) 


45 


4.494 (2) 


2S1/2 


-243 


55937 





04627 


-0.00856 


-0.00034 


0.00006 


± 0.00002 


-243 


52194 


(3) 






2P1/2 


-240 


16111 





21728 


-0.00023 


-0.00011 


-0.00010 


± 0.00003 


-239 


94427 


(3) 






2P3/2 


-233 


81565 





08940 


-0.00690 


0.00066 


0.00002 


± 0.00004 


-233 


73247 


(5) 


46 


4.532 (3) 


2Si/ 2 


-255 


32827 





05205 


-0.00873 


-0.00036 


0.00006 


± 0.00002 


-255 


28526 


(3) 






2P1/2 


-251 


83223 





23411 


-0.00025 


-0.00011 


-0.00010 


± 0.00003 


-251 


59857 


(4) 






2p 3 /2 


-244 


85971 





09595 


-0.00710 


0.00073 


0.00002 


± 0.00005 


-244 


77012 


(5) 


47 


4.544 (4) 


2Si/ 2 


-267 


40564 (1) 





05814 


-0.00891 


-0.00038 


0.00006 


± 0.00003 


-267 


35673 


(3) 






2P1/2 


-263 


81030 





25179 


-0.00027 


-0.00011 


-0.00010 


± 0.00004 


-263 


55899 


(4) 






2P3/2 


-256 


16366 





10276 


-0.00731 


0.00080 


0.00002 


± 0.00005 


-256 


06738 


(5) 


48 


4.614(2) 


2*1/2 


-279 


79438 (1) 





06459 


-0.00909 


-0.00040 


0.00006 


± 0.00003 


-279 


73921 


(3) 






2j5l/2 


-276 


09840 





27037 


-0.00030 


-0.00011 


-0.00010 


± 0.00004 


-275 


82853 


(4) 






2P3/2 


-267 


72796 





10984 


-0.00752 


0.00089 


0.00002 


± 0.00006 


-267 


62473 


(6) 


49 


4.617(2) 


2Si/ 2 


-292 


49801 (1) 





07139 


-0.00927 


-0.00042 


0.00006 


± 0.00003 


-292 


43624 


(3) 






2P1/2 


-288 


69965 





28988 


-0.00033 


-0.00011 


-0.00010 


± 0.00004 


-288 


41032 


(4) 






2P3/2 


-279 


55310 





11720 


-0.00774 


0.00097 


0.00002 


± 0.00006 


-279 


44265 


(6) 


50 


4.654(1) 


2Si/ 2 


-305 


51949 (1) 





07853 


-0.00946 


-0.00045 


0.00007 


± 0.00003 


-305 


45080 


(3) 






2pi/ 2 


-301 


61727 





31030 


-0.00038 


-0.00010 


-0.00011 


± 0.00005 


-301 


30756 


(5) 






2P3/2 


-291 


63957 





12483 


-0.00797 


0.00107 


0.00002 


± 0.00006 


-291 


52162 


(7) 


51 


4.680 (2) 


2Si/ 2 


-318 


86248 (2) 





08610 


-0.00965 


-0.00047 


0.00007 


± 0.00003 


-318 


78643 


(4) 






2pi/ 2 


-314 


85471 





33175 


-0.00043 


-0.00010 


-0.00011 


± 0.00005 


-314 


52359 


(5) 






2P3/2 


-303 


98789 





13275 


-0.00821 


0.00117 


0.00002 


± 0.00007 


-303 


86216 


(7) 


52 


4.743 (3) 


2Si/ 2 


-332 


53029 (2) 





09404 


-0.00985 


-0.00049 


0.00007 


± 0.00003 


-332 


44652 


(4) 






2^1/2 


-328 


41543 





35420 


-0.00048 


-0.00010 


-0.00011 


± 0.00005 


-328 


06192 


(5) 






2P3/2 


-316 


59857 





14094 


-0.00845 


0.00128 


0.00001 


± 0.00007 


-316 


46479 


(7) 


53 


4.750 (4) 


2S1/2 


-346 


52697 (3) 





10232 


-0.01005 


-0.00051 


0.00007 


± 0.00003 


-346 


43514 


(5) 






2P1/2 


-342 


30302 





37764 


-0.00055 


-0.00009 


-0.00011 


± 0.00005 


-341 


92614 


(6) 






2P3/2 


-329 


47213 





14942 


-0.00871 


0.00139 


0.00001 


± 0.00008 


-329 


33002 


(8) 


54 


4.787(5) 


2Si/ 2 


-360 


85598 (4) 





11103 


-0.01026 


-0.00054 


0.00007 


± 0.00004 


-360 


75567 


(5) 






2pi/ 2 


-356 


52130 





40217 


-0.00062 


-0.00009 


-0.00012 


± 0.00006 


-356 


11995 


(6) 






2P3/2 


-342 


60914 





15818 


-0.00897 


0.00152 


0.00001 


± 0.00008 


-342 


45840 


(8) 


55 


4.804(5) 


2sr/2 


-375 


52155 (4) 





12016 


-0.01047 


-0.00056 


0.00008 


± 0.00004 


-375 


41234 


(6) 






2P1/2 


-371 


07417 





42782 


-0.00070 


-0.00008 


-0.00012 


± 0.00006 


-370 


64725 


(6) 






2P3/2 


-356 


01014 





16723 


-0.00924 


0.00165 


0.00001 


± 0.00009 


-355 


85050 


(9) 


56 


4.838 (5) 


2Si/ 2 


-390 


52754 (5) 





12971 


-0.01069 


-0.00058 


0.00008 


± 0.00004 


-390 


40902 


(6) 






2P1/2 


-385 


96566 





45460 


-0.00079 


-0.00008 


-0.00012 


± 0.00007 


-385 


51205 


(7) 






2P3/2 


-369 


67571 





17656 


-0.00952 


0.00179 


0.00000 


± 0.00009 


-369 


50688 


(10) 


57 


4.855 (5) 


2si/ 2 


-405 


87846 (6) 





13972 


-0.01091 


-0.00061 


0.00008 


± 0.00004 


-405 


75018 


(7) 






2^1/2 


-401 


19999 





48257 


-0.00088 


-0.00007 


-0.00013 


± 0.00007 


-400 


71850 


(7) 






2^3/2 


-383 


60643 





18618 


-0.00981 


0.00194 


0.00000 


± 0.00010 


-383 


42812 


(10) 


58 


4.877(2) 


2Si/ 2 


-421 


57854 (4) 





15015 


-0.01114 


-0.00063 


0.00008 


± 0.00004 


-421 


44008 


(6) 




2^1/2 


-416 


78149 





51173 


-0.00099 


-0.00006 


-0.00013 


± 0.00007 


-416 


27094 


(7) 






2P3/2 


-397 


80290 





19609 


-0.01011 


0.00209 


0.00000 


±0.00011 


-397 


61482 


(11) 


59 


4.892 (5) 


2Si/ 2 


-437 


63243 (7) 





16103 


-0.01138 


-0.00065 


0.00008 


± 0.00005 


-437 


48335 


(9) 






2P1/2 


-432 


71467 





54213 


-0.00110 


-0.00005 


-0.00013 


± 0.00008 


-432 


17382 


(8) 






2P3/2 


-412 


26571 





20629 


-0.01042 


0.00226 


0.00000 


±0.00011 


-412 


06758 


(11) 


60 


4.912(2) 


2Si/ 2 


-454 


04479 (6) 





17240 


-0.01162 


-0.00068 


0.00009 


± 0.00005 


-453 


88459 


(8) 






2pi/ 2 


-449 


00422 





57384 


-0.00122 


-0.00004 


-0.00014 


± 0.00008 


-448 


43178 


(8) 






z P3/2 


-426 


99551 





21678 


-0.01073 


0.00244 


— n nnnm 


4- o nnm 9 


-426 


78703 




61 


4.962 (50) 


2Si/ 2 


-470 


82019 (69) 





18416 


-0.01187 


-0.00070 


0.00009 


± 0.00005 


-470 


64850 


(70) 






2P1/2 


-465 


65499 (3) 





60677 


-0.00136 


-0.00003 


-0.00014 


± 0.00009 


-465 


04975 


(9) 






2P3/2 


-441 


99292 





22756 


-0.01106 


0.00262 


-0.00001 


± 0.00013 


-441 


77382 


(13) 


62 


5.084(6) 


2Si/ 2 


-487 


96303 (11) 





19662 


-0.01212 


-0.00072 


0.00009 


± 0.00005 


-487 


77916 


(12) 






2pi/ 2 


-482 


67204 





64125 


-0.00150 


-0.00002 


-0.00014 


± 0.00009 


-482 


03246 


(9) 






2P3/2 


-457 


25862 





23862 


-0.01141 


0.00282 


-0.00001 


± 0.00013 


-457 


02860 


(13) 



TABLE VI: Electronic-structure contributions (continued) 



z 


<rV' a 


State 


Dirac 


1 -photon 


2ph.[MBPT] 2ph.[QED] 


3-photon 




h.o. 


Total 




63 


5.113(4) 


2si/ 2 


-505 


48081 (10) 





20935 


-0.01238 


-0 


00075 





00009 


± 





00005 


-505 


28450 


(11) 






2Pl/2 


-500 


06064 





67693 


-0.00164 


-0 


00001 


-0 


00014 


± 





00010 


-499 


38550 


(10) 






2P3/2 


-472 


79322 





24997 


-0.01176 





00303 


-0 


00002 


± 





00014 


-472 


55200 


(14) 


64 


5.162 (2) 


2S1/2 


-523 


37749 (9) 





22275 


-0.01265 


-0 


00077 





00009 


± 





00006 


-523 


16807 


(11) 






2P1/2 


-517 


82618 





71421 


-0.00180 





00000 


-0 


00015 


± 





00010 


-517 


11392 


(10) 






2P3/2 


-488 


59744 





26161 


-0.01212 





00325 


-0 


00002 


± 





00015 


-488 


34473 


(15) 


65 


5.060(150) 


2S1/2 


-541 


6620 (32) 





23657 


-0.01293 


-0 


00079 





00010 


± 





00006 


-541 


4391 (32) 






2P1/2 


-535 


97441 (16) 





75293 


-0.00197 





00001 


-0 


00015 


± 





00011 


-535 


22358 


(19) 






2P3/2 


-504 


67192 





27353 


-0.01249 





00348 


-0 


00002 


± 





00015 


-504 


40743 


(15) 


66 


5.221 (2) 


2Si/ 2 


-560 


33217(13) 





25100 


-0.01321 


-0 


00081 





00010 


± 





00006 


-560 


09510 


(14) 






2P1/2 


-554 


51094 





79318 


-0.00215 





00003 


-0 


00016 


± 





00012 


-553 


72003 


(12) 






2P3/2 


-521 


01747 





28573 


-0.01288 





00372 


-0 


00003 


± 





00016 


-520 


74092 


(16) 


67 


5.202 (31) 


2S1/2 


-579 


40308 (83) 





26596 


-0.01350 


-0 


00083 





00010 


± 





00006 


-579 


15136 


(83) 






2P1/2 


-573 


44216 (4) 





83506 


-0.00234 





00004 


-0 


00016 


± 





00012 


-572 


60956 


(13) 






2P3/2 


-537 


63469 





29821 


-0.01328 





00398 


-0 


00003 


± 





00017 


-537 


34581 


(17) 


68 


5.250 (3) 


2*1/2 


-598 


87578 (17) 





28145 


-0.01380 


-0 


00085 





00010 


± 





00007 


-598 


60888 


(18) 






2P1/2 


-592 


77424 





87853 


-0.00254 





00006 


-0 


00016 


± 





00013 


-591 


89835 


(13) 






2P3/2 


-554 


52439 





31097 


-0.01369 





00425 


-0 


00004 


± 





00018 


-554 


22290 


(18) 


69 


5.226 (4) 


2Si/ 2 


-618 


76050 (20) 





29754 


-0.01411 


-0 


00087 





00011 


± 





00007 


-618 


47784 


(21) 






2P1/2 


-612 


51384(1) 





92373 


-0.00274 





00007 


-0 


00017 


± 





00014 


-611 


59295 


(14) 






2P3/2 


-571 


68728 





32401 


-0.01412 





00453 


-0 


00004 


± 





00019 


-571 


37290 


(19) 


70 


5.311(6) 


2S1/2 


-639 


05870 (28) 





31424 


-0.01443 


-0 


00089 





00011 


± 





00007 


-638 


75966 


(29) 






2P1/2 


-632 


66766 (2) 





97068 


-0.00296 





00009 


-0 


00017 


± 





00014 


-631 


70003 


(14) 






2P3/2 


-589 


12417 





33731 


-0.01456 





00482 


-0 


00005 


± 





00019 


-588 


79664 


(19) 


71 


5.370 (30) 


2S1/2 


-659 


7814(12) 





33151 


-0.01476 


-0 


00090 





00011 


± 





00008 


-659 


4655 (12) 






2P1/2 


-653 


24294 (7) 


1 


01940 


-0.00319 





00011 


-0 


00017 


± 





00015 


-652 


22680 


(17) 






2P3/2 


-606 


83581 





35088 


-0.01501 





00514 


-0 


00005 


± 





00020 


-606 


49485 


(20) 


72 


5.342 (2) 


2Si/ 2 


-680 


93893 (26) 





34931 


-0.01509 


-0 


00092 





00011 


± 





00008 


-680 


60551 


(28) 






2P1/2 


-674 


24721 (2) 


1 


06994 


-0.00344 





00013 


-0 


00018 


± 





00016 


-673 


18076 


(16) 






2P3/2 


-624 


82296 





36472 


-0.01547 





00546 


-0 


00006 


± 





00021 


-624 


46831 


(21) 


73 


5.351(3) 


2S1/2 


-702 


53427 (32) 





36780 


-0.01544 


-0 


00093 





00012 


± 





00008 


-702 


18272 


(33) 






2P1/2 


-695 


68798 (2) 


1 


12246 


-0.00369 





00014 


-0 


00018 


± 





00017 


-694 


56925 


(17) 






2P3/2 


-643 


08648 





37881 


-0.01595 





00580 


-0 


00006 


± 





00022 


-642 


71788 


(22) 


74 


5.367(2) 


2Si/ 2 


-724 


57660 (32) 





38692 


-0.01579 


-0 


00094 





00012 


± 





00008 


-724 


20630 


(33) 






2pi/ 2 


-717 


57343 (2) 


1 


17695 


-0.00395 





00016 


-0 


00019 


± 





00017 


-716 


40046 


(18) 






2P3/2 


-661 


62718 





39316 


-0.01645 





00616 


-0 


00007 


± 





00023 


-661 


24438 


(23) 


75 


5.339(13) 


2sr/2 


-747 


07707 (85) 





40664 


-0.01616 


-0 


00095 





00012 


± 





00009 


-746 


68742 


(85) 






2P1/2 


-739 


91201 (6) 


1 


23347 


-0.00423 





00018 


-0 


00019 


± 





00018 


-738 


68278 


(19) 






2P3/2 


-680 


44589 





40775 


-0.01696 





00653 


-0 


00007 


± 





00024 


-680 


04863 


(24) 


76 


5.413(1) 


2S1/2 


-770 


03598 (40) 





42704 


-0.01654 


-0 


00096 





00013 


± 





00009 


-769 


62631 


(41) 






2P1/2 


-762 


71201 (3) 


1 


29208 


-0.00452 





00020 


-0 


00020 


± 





00019 


-761 


42444 


(19) 






2P3/2 


-699 


54352 





42259 


-0.01748 





00692 


-0 


00008 


± 





00025 


-699 


13157 


(25) 


77 


5.402 (106) 


2Si/ 2 


-793 


4734 (80) 





44808 


-0.01693 


-0 


00097 





00013 


± 





00009 


-793 


0431 (80) 






2P1/2 


-785 


98315 (59) 


1 


35290 


-0.00482 





00022 


-0 


00020 


± 





00020 


-784 


63506 


(62) 






2P3/2 


-718 


92087 





43767 


-0.01802 





00733 


-0 


00009 


± 





00026 


-718 


49398 


(26) 


78 


5.428 (3) 


2Si/ 2 


-817 


39139 (55) 





46992 


-0.01733 


-0 


00097 





00013 


± 





00010 


-816 


93964 


(56) 






2P1/2 


-809 


73458 (4) 


1 


41607 


-0.00513 





00024 


-0 


00021 


± 





00021 


-808 


32362 


(21) 






2P3/2 


-738 


57887 





45297 


-0.01858 





00775 


-0 


00009 


± 





00027 


-738 


13682 


(27) 


79 


5.436 (4) 


2Si/ 2 


-841 


80344 (67) 





49238 


-0.01774 


-0 


00097 





00014 


± 





00010 


-841 


32964 


(67) 






2P1/2 


-833 


97640 (5) 


1 


48152 


-0.00546 





00027 


-0 


00021 


± 





00022 


-832 


50029 


(23) 






2P3/2 


-758 


51842 





46850 


-0.01916 





00820 


-0 


00010 


± 





00028 


-758 


06098 


(28) 


80 


5.463 (2) 


2Si/ 2 


-866 


71677 (67) 





51557 


-0.01817 


-0 


00097 





00014 


± 





00010 


-866 


22020 


(67) 






2P1/2 


-858 


71879 (5) 


1 


54942 


-0.00580 





00029 


-0 


00022 


± 





00023 


-857 


17510 


(24) 






2P3/2 


-778 


74045 





48423 


-0.01975 





00866 


-0 


00011 


± 





00029 


-778 


26741 


(29) 


81 


5.476 (3) 


2Si/ 2 


-892 


14506 (79) 





53945 


-0.01861 


-0 


00096 





00014 


± 





00011 


-891 


62504 


(79) 






2pi/ 2 


-883 


97272 (6) 


1 


61981 


-0.00615 





00031 


-0 


00023 


± 





00024 


-882 


35898 


(25) 






2P3/2 


-799 


24589 





50018 


-0.02035 





00915 


-0 


00011 


± 





00030 


-798 


75703 


(30) 


82 


5.501(1) 


2Si/ 2 


-918 


09688 (83) 





56397 


-0.01906 


-0 


00096 





00015 


± 





00011 


-917 


55279 


(84) 






2P1/2 


-909 


74938 (7) 


1 


69271 


-0.00652 





00033 


-0 


00023 


± 





00025 


-908 


06309 


(26) 






2P3/2 


-820 


03570 





51632 


-0.02098 





00965 


-0 


00012 


± 





00031 


-819 


53083 


(31) 



20 



TABLE VI: Electronic-structure contributions (continued) 



z 


<rV /2 


State 


Dirac 


1 -photon 


2ph.[MBPT] 2ph.[QED] 


3-photon 




h.o. 


Total 


83 


5.521(3) 


2Si/ 2 


-944 


.5857(10) 





.58922 


-0 


.01953 


-0 


.00095 





.00015 


± 





.00012 


-944 


.0168 (10) 






2P1/2 


-936 


.06078 (9) 


1 


.76830 


-0 


.00690 





.00035 


-0 


.00024 


± 





.00026 


-934 


.29927 (28) 






2P3/2 


-841 


.11086 





.53265 


-0 


.02162 





.01017 


-0 


.00013 


± 





.00032 


-840 


.58979 (32) 


84 


5.526(13) 


2S1/2 


-971 


.6251 (22) 





.61524 


-0 


.02002 


-0 


.00093 
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± 





.00012 


-971 


.0307 (22) 






2P1/2 


-962 


.91947 (20) 


1 


.84671 


-0 


.00730 





.00037 


-0 


.00024 


± 





.00028 


-961 


.07993 (34) 






2P3/2 


-862 


.47235 





.54915 


-0 


.02229 





.01072 


-0 


.00014 


± 





.00034 


-861 


.93490 (34) 


85 


5.539 (55) 


2S1/2 


-999 


.2252 (94) 





.64207 


-0 


.02052 


-0 


.00092 





.00016 


± 





.00012 


-998 


.6044 (94) 






2P1/2 


-990 


.33823 (89) 


1 


.92808 


-0 


.00771 





.00039 


-0 


.00025 


± 





.00029 


-988 


.41772 (93) 






2P3/2 


-884 


.12119 





.56583 


-0 


.02297 





.01129 


-0 


.00014 


± 





.00035 


-883 


.56719 (35) 


86 


5.655(16) 


2Si/ 2 


-1027. 


.3822 (33) 





.66963 


-0 


.02104 


-0 


.00090 





.00016 


± 





.00013 


-1026 


.7343 (33) 






2P1/2 


-1018 


.32939 (33) 


2 


.01237 


-0 


.00814 





.00041 


-0 


.00026 


± 





.00030 


-1016 


.32500 (44) 






2f>3/2 


-906 


.05854 





.58266 


-0 


.02367 





.01188 


-0 


.00015 


± 





.00036 


-905 


.48783 (36) 


87 


5.658(13) 


2si/ 2 


-1056 


.1458 (31) 





.69798 


-0 


.02157 


-0 


.00087 





.00017 


± 


0. 


.00013 


-1055 


.4701 (31) 






2P1/2 


-1046. 


.91008 (31) 


2 


.09987 


-0 


.00858 





.00043 


-0 


,00026 


± 


0. 


00031 


-1044 


.81862 (44) 






2^3/2 


-928 


.28519 





.59963 


-0 


.02439 





.01249 


-0 


.00016 


± 


0. 


.00037 


-927 


.69762 (37) 


88 


5.684 (26) 


2Sl/2 


-1085. 


.5082 (62) 





.72710 


-0 


.02212 


-0 


.00084 





.00017 


± 


0. 


.00014 


-1084 


.8039 (62) 






2pi/ 2 


-1076 


.09350 (65) 


2 


.19058 


-0 


.00904 





.00045 


-0 


.00027 


± 


0. 


.00033 


-1073 


.91178 (73) 






2P3/2 


-950 


.80236 





.61674 


-0 


.02513 





.01313 


-0 


.00017 


± 


0. 


.00038 


-950 


.19780 (38) 


89 


5.670 (57) 


2Si/ 2 


-1115. 


.498 (14) 





.75687 


-0 


.02269 


-0 


.00081 





.00017 


± 


0. 


.00014 


-1114 


.764 (14) 






2P1/2 


-1105. 


.8966(16) 


2 


.28455 


-0 


.00951 





.00047 


-0 


,00028 


± 


0. 


.00034 


-1103 


.6214 (16) 






2P3/2 


-973 


.61105 





.63395 


-0 


.02590 





.01379 


-0 


.00018 


± 


0. 


.00039 


-972 


,98938 (39) 


90 


5.710(50) 


2Si/ 2 


-1146. 


.109(14) 





.78750 


-0 


.02328 


-0 


.00077 





,00018 


± 


0. 


.00015 


-1145 


.346 (14) 
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2^3/2 


-996 


.71252 





.65128 


-0 


.02668 





.01448 


-0 


.00019 


± 


0. 


.00041 


-996 


.07363 (41) 


91 


5.700 (57) 


2Si/ 2 


-1177. 


.385(18) 





.81929 


-0 


.02385 


-0 


.00072 





.00018 


± 


0. 


.00015 


-1176 


.590(18) 






2pi/ 2 


-1167. 


.4260 (20) 


2 


.48333 


-0 


.01052 





.00051 


-0 


.00029 


± 





.00037 


-1164 


.9529(21) 






2P3/2 


-1020. 


.10774 





.66872 


-0 


.02749 





.01519 


-0 


.00020 


± 


0. 


.00042 


-1019 


.45151 (42) 


92 


5.851 (7) 


2si/ 2 


-1209. 


.2997(37) 





.87286 


-0 


.02466 


-0 


.00067 





.00019 


± 


0. 


.00016 


-1208 


.4520 (37) 






2P1/2 


-1199. 


.18416 (42) 


2. 


.58825 


-0 


.01104 





.00053 


-0 


.00030 


± 


0. 


.00038 


-1196 


.60672 (57) 






2f>3/2 


-1043. 


.79820 





.68624 


-0 


.02831 





.01593 


-0 


.00021 


± 


0. 


.00043 


-1043 


.12456 (43) 



